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PHASE -CONJUGATED  FLUORESCENCE 


Henk  F.  Arnoldus 
Department  of  Physics 
Mendel  Hall 
Villanova  University 
Villanova,  Pennsylvania  19085 

and 

Thomas  F.  George 

Departments  of  Chemistry  and  Physics  &  Astronomy 
239  Fronczak  Hall 

State  University  of  New  York  at  Buffalo 
Buffalo,  New  York  14260 

Fluorescent  emission  by  an  atom  near  a  phase  conjugator  (PC)  based  on 
four-wave  mixing  is  studied  from  first  principles.  The  Maxwell-Heisenberg 
equations  are  solved  for  the  radiation  field,  and  with  an  asymptotic  expansion 
an  expression  is  derived  for  the  field  in  the  far  zone.  The  total  emitted 
power  which  can  be  measured  by  a  detector  in  this  region  is  evaluated,  and  it 
is  found  that  this  power  acquires  three  distinct  contributions.  First,  there 
are  photons  which  are  emitted  by  the  atom  directly  towards  the  detector,  and 
without  any  interaction  with  the  medium.  Second,  there  are  photons  that  first 
travel  towards  the  surface  of  the  PC,  and  they  have  a  certain  probability  of 
being  reflected  in  the  specular  direction  and  towards  the  detector.  The  third 
kind  of  radiation  consists  of  phase - conj ugated  photons,  which  are  emitted 
independently  of  the  previous  ones.  It  is  shown  that  the  first  two  processes 
are  a  result  of  simple  atomic  spontaneous  decay,  but  that  the  emmission  of  a 
phase - conj ugated  fluorescent  photon  involves  a  three-photon  process.  The 
latter  process  has  a  probability  proportional  to  the  population  of  the  atomic 
ground  state.  It  is  pointed  out  that  an  atom  in  its  ground  state  polarizes 
the  nonlinear  medium  of  the  PC,  which  subsequently  can  emit  spontaneously  two 
photons.  An  absorption-emission- absorption  process  by  the  atom  then  produces 
a  fluorescent  photon,  together  with  a  spontaneous  excitation  of  the  atom. 


PACS:  42 .65 .H,  32.80,  52.50.D 
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I.  INTRODUCTION 

Wavefront  inversion  of  optical  radiation  was  observed  for  the  first  time 
by  Zel'dovich  et  al ^  in  stimulated  Brillouin  scattering.  Mainly  due  to  the 

fact  that  wavefront  inversion,  or  phase  conjugation,  can  be  utilized  to 

2  3 

correct  wavefront  distortions  in  optical  amplifiers,  ’  this  technique  has 
since  attracted  a  considerable  amount  of  attention.  In  1977  it  was  proposed 
by  Hellwarth  and  Yariv  and  Pepper  to  construct  a  phase  conjugator  (PC)  which 
operates  through  four-wave  mixing  (FWM)  in  a  nonlinear  medium,  and  almost 
simultaneously  this  was  realized  experimentally  in  a  liquid  CS2  cell^’^  and  in 

g 

a  lithium  formate  crystal.  Alternative  methods  of  phase  conjugation  include 
9-11  1213 

FWM  in  gases  and  thin  films,  ’  or  more  exotic  media  like  microparticles 

-  •  .  14  15 

suspended  in  a  liquid*  or  organic  dye  molecules  in  a  solid  matrix.  Also, 

the  change  in  reflectivity  of  a  thin  metal  layer  under  high  irradiance  has 

been  considered  as  a  candidate  for  the  generation  of  phase-conjugated 

light.  ’  More  recently,  phase  conjugation  in  photorefractive  crystals, 

1  Q 

like  BaTiO^,  has  gained  popularity,  which  can  be  attributed  to  the 

possibility  of  self-pumping  of  the  crystal  by  the  incident  beam  whose  phase- 

conjugated  image  is  sought.  ^  ^  On  the  theoretical  side,  much  effort  has 

been  devoted  to  the  description  of  (idealized)  phase-conjugated  radiation,  and 
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to  the  study  of  the  applicability  of  a  PC  for  wavefront  correction. 

The  possibility  of  wavefront  distortion  correction  after  phase 
conjugation  relies  on  the  equivalence  between  phase  conjugation  and  time 
reversal,  as  explained  in  various  reviews.  ’  Taking  the  complex  conjugate 
of  the  spatial  part  (the  phase)  of  a  field  is  mathematically  equivalent  to 
replacing  t  by  -t.  Many  of  the  properties  of  phase-conjugated  radiation  can 

be  understood  most  easily  with  time- reversal  arguments.  For  instance,  when  a 

plane  wave  is  incident  on  a  phase-conjugating  crystal  (under  a  certain  angle), 


chen  the  time -reversed  replica  of  this  beam  must  again  be  a  plane  wave,  and  it 


propagates  in  the  direction  opposite  to  the  incident  wave.  This  is  in 
contrast  to  specular  reflection  at  a  linear  medium  (dielectric  or  metal). 
Also,  when  a  point  source  of  radiation,  like  a  fluorescing  atom,  is  close  to 
the  surface  of  a  PC,  then  the  incident  field  is  a  diverging  spherical  wave. 
According  to  the  time-reversal  argument,  the  reflected  (phase-conjugated)  wave 
must  be  a  converging  spherical  wave  which  is  focused  exactly  on  the  point 
source.  This  reflected  field  at  the  position  of  the  source  can  affect  the 

dynamical  evolution  of  the  radiatior  considerably,  and  thereby  its  mechanism 
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of  emission,  as  was  realized  for  the  first  time  by  Agarwal. 

Recently,  Sochove''*  has  redone  the  calculations  of  Ref.  30,  with  an 

approach  where  the  radiation  field  is  quantized  explicitly.  In  both 

references  the  Einstein  coefficient  for  spontaneous  decay  of  an  oscillating 

dipole  was  obtained,  and  it  appeared  indeed  that  there  are  some  remarkable 

differences  as  compared  to  fluorescent  emission  in  empty  space  or  near  a  metal 

surface.  Most  notably,  it  was  found  that  the  decay  race  is  independent  of  the 

distance  between  the  dipole  and  the  PC,  as  could  be  expected  from  the 

equivalence  of  phase  conjugation  and  time  reversal.  A  reflected  photon 

returns  to  the  atoms,  no  matter  how  far  away  this  atom  is.  Of  course,  there 

are  some  limitations  on  this  result  due  to  causality  requirements  and 

retardation  effects,  as  recognized  by  the  authors.  Recently,  Cook  and 

Milonni  have  shown  that  a  sample  of  many  two- level  atoms  in  their  ground 

states  and  near  a  PC  is  unstable.  We  shall  show  that  this  result  also  holds 

for  a  single  atom,  and  explain  the  underlying  physical  mechanism.  Hendricks 
3  3 

and  NienhuisJ  have  also  studied  the  spontaneous  decay  of  a  two-level  atom 
near  a  PC.  Their  results  for  the  decay  rate  are  consistent  with  our  findings 
for  the  emission  rate. 


j 


4 


We  consider  an  atomic  dipole  p(t)  with  arbitrary  time  dependence  and 
evaluate  the  radiation  field  in  the  far  zone,  as  it  can  be  observed 
experimentally  by  a  photodetector.  The  time  dependence  of  the  Heisenberg 
operator  /i(t)  is  brought  about  by  its  interaction  with  all  components  of  the 
electric  field.  We  were  able  to  avoid  an  explicit  quantization  in  terms  of 
plane-wave  modes  of  the  radiation  field,  and  the  general  results  do  not  depend 
in  form  on  the  specific  properties  of  the  dipole  moment  n(t) .  We  work  out  the 
case  of  a  two- level  atom,  and  it  appears  that  some  quite  remarkable  features 
in  the  behavior  of  an  atom  near  a  PC  can  be  predicted  without  knowledge  of  the 
details  of  the  temporal  evolution  of  the  atomic  density  operator. 

II .  VARIOUS  FIELD  COMPONENTS 

A  nonlinear  medium  occupies  a  part  of  the  region  z  <  0,  and  its  surface 
is  the  plane  z  -  0.  Two  counterpropagating  laser  beams  with  frequency  u  >  0 
pump  the  medium,  and  the  third-order  nonlinear  susceptibility  is  responsible 
for  a  four-wave  mixing  process  between  the  two  pump  fields,  an  incident  field, 
and  a  generated  phase  -  conj  ugated  wave.  We  shall  allow  the  medium  to  have  a 
non-unity  dielectric  constant,  which  implies  that  we  have  to  take  into  account 
the  specular  reflection  of  an  incident  wave.  In  this  fashion  we  can  keep 
track  of  the  differences  between  linear  and  nonlinear  effects,  since  they  are 
both  unified  in  a  single  formalism.  A  dipole  /i(t)  is  located  at  r  -  hez>  h  > 
0,  on  the  positive  z-axis. 

We  are  interested  in  the  electric  field  operator  E(r,t)  in  the  region 
z  >  0  only,  since  that  is  the  place  where  we  can  put  a  detector  and  measure 
the  radiation.  In  general,  E(r,t)  acquires  three  distinct  contributions,  and 


we  can  write 
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E(r,t)  =  Ef(r,t)  +  EQ(r.t)  +  E^r.t)  (2.1) 

for  the  field  in  the  Heisenberg  picture.  Here,  is  the  free  field  which 

includes  the  vacuum  field  and  possible  external  fields.  The  term  E  is  the 

~P 

particular  solution  and  is  equal  to  the  radiation  field  of  an  atomic  dipole  in 
empty  space.  A  part  of  this  field  is  incident  upon  the  PC,  which  generates  a 
phase-conjugated  signal  and  a  reflected  specular  wave,  both  of  which  are 
included  in  E.  (h  —  homogeneous  solution) . 

The  key  to  the  evaluation  of  E(r,t)  is  the  notion  that  the  Heisenberg 
equations  for  the  time  evolution  of  the  e  lectromagnetic  -  field  operators  are 
identical  in  form  to  the  classical  Maxwell's  equations  .  Just  as  in 

classical  electrodynamics,  we  cannot  consider  the  field  in  the  region  z  >  0 
only,  but  we  have  to  take  into  account  the  solution  in  z  <  0  as  well.  At  the 
boundary  the  two  solutions  must  match  in  the  usual  way.  The  free  field  E^  is 
assumed  to  obey  Maxwell's  equations  separately.  For  parametric  FWM ,  these 
equations  are  linear  in  the  fields  (in  the  undepleted-pump  approximation),  and 
therefore  E-Ep  must  also  be  a  solution  of  Maxwell's  equations.  We  call  this 
the  source  field  and  write  it  as 

5a(£.=)  “  E  (r,t)  +  E.n(r.t)  ,  (2.2) 

for  z  >  0.  This  component  of  the  radiation  field  is  due  to  the  presence  of 
the  dipole  and  its  interaction  with  the  medium.  We  shall  focus  our  attention 
on  Eg  only,  since  this  is  the  radiative  part  which  can  be  measured  by  a 
detector  in  the  far  zone.  The  homogeneous  component  has  two  distinct 


contributions , 
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Eh(r,c)  -  Er(r,t)  +  E  £,t) 


(2.3) 


in  terms  of  a  specular  (r)  and  a  phase-conjugated  (pc)  part. 

It  is  convenient  to  adopt  a  Fourier  transform  of  E(r,t)  according  to 


*->  -  r. 


icjt 

dt  e  E(r , t) 


(2.4) 


The  electric-field  operator  is  Hermitian,  which  translates  into  the  Fourier 
domain  as 


E(r ,  t) 1  -  E(r,t)  — •  E  ( r  ,  _j  ;  -  E(r,-w) 


(2.5) 


A  useful  concept  is  the  positive- frequency  part  of  the  field,  which  is  defined 


E(r , t) (+)  -  J-  |  dw 

J  n 


'iwt  p/  \ 
e  E(r,w) 


(2.6) 


in  terms  of  E(r,u).  The  negative  -  frequency  part  follows  from 


E(r,t)(')  -  [E(r,c)(+)]t 


(2.7) 


and  the  total  field  assumes  the  form 


E(r,t)  -  E(r,t)(+)  +  E(r,t)('} 


(2.8) 
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Similar  notations  will  be  used  for  other  time  -  dependent  quantities,  and  most 
notably  for  the  dipole  operator  -(.t).  which  is  also  Hermitian.  From  Eq.  (2.5) 
it  follows  that  we  know  the  field  E(r,o)  as  soon  as  its  value  for  w  >  0  only 
is  given.  Zauivalentl  y ,  Eqs  .  (2.7)  and  (2.8)  show  that  the  entire  field  is 
determined  by  just  E''  +  ')  or  E^  \ 


III.  ANGULAR  SPECTRUM  OF  PLANE  WAVES 


In  this  section  we  derive  an  explicit 
for  a  dipole  in  empty  space,  which  is  most 
homogeneous  contribution  i'Sec.  V)  .  For  w  > 
electric  field  of  a  dipole  at  r  -  h  assumes 


form  of  the  particular  solution  E^ 
suitable  for  the  evaluation  of  the 
0,  the  Fourier  transform  of  the 
the  form 


47T  i 


(k2a(w)  +  (p(w).V)V)  ~ 


ik I r-h | 
1  r-h| 


0.1) 


with  k  -  w/c  ,  and  where  only  a(w)  is  a  quantum  operator.  After  carrying  out 
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the  differentiations,  this  expression  reduces  to  the  more  familiar  form.  We 

36 

shall  use  Wevl's  representation  of  the  Green's  function" 


e^lr-hl 


63  -  e  iax+i^y+i-y  |  z  - h  | 

"V 

J  -CO 


k  >  0 


(3.2) 


The  parameter  7  is  defined  by 


>  7 

-  3" 


7  “ 


(3.3) 
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and  we  cake  the  form  for  which  the  argument  of  the  square  root  is  positive. 

Substituting  expression  (3.2)  into  Eq.  (3.1)  and  carrying  out  the 
differentiations  yields 


E  (r,w) 


^  M(w)j6 (r-h) 


8jt  t  J  -co 
o 


j>r 


1  iax+i^y+i7 | z-h j 


dfl  A  e 


x  Ik  u(w)  -  .Q^^(w)  +  +  -rsgn(z-h)M,(w) 


x  [ae  +  0e  +  7Sgn(z-h)e  ] ) 


(3.4) 


Here,  m(w)j_  i-s  the  perpendicular  component  of  n(u)  with  respect  to  the  plane 
z  "  0.  The  <5 -function  at  r  -  h  appears  due  to  the  second  derivative  of  |z-h| 
with  respect  to  z.  For  two  given  values  of  the  integration  variables  a  and  f), 
we  define  the  vector 


Kii  “  ae  +  /3e 
~x  -y 


(3.5) 


where  the  subscript  ]j  indicates  that  this  vector  is  parallel  to  the  z  —  0 
plane.  In  terms  of  7,  we  then  introduce  the  two  complementary  vectors 


5  *  h  * 


(3.6) 


5’  ‘  h  ■  nz  ■ 


(3.7) 
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which  are  each  other's  mirror  image.  Notice  that  Ki|  is  real,  but  that  K  and 
K'  can  each  have  an  imaginary  z -component.  For  the  square  of  the  magnitude  of 
the  various  wave  vectors,  we  find 


1  0  0  0  0  1 

K-  -  (K* )  -  k~  -  a"  +  d  +  7‘ 


(«/c)2 


(3.8) 


Furthermore,  we  write  symbolically 


2 

r 

d  K;|  (  .  ■ 

.  • )  -  da 

J  -CO 

dB  (...)  . 


(3.9) 


With  the  notations  from  above,  E  (r,w)  can  be  cast  in  the  form 

-p  - 


E  (r.y)  -  — \ —  [  d2K„  1  eL-*(S’-) 
-P  ~  8?r2e  J  ~:l  1 


x  |k*>(u)  -  i^(«)*K)K),  for  z  >  h 


(3.10) 


and  the  expression  for  0  <  z  <  h  follows  from  Eq  .  (3.10)  after  the 

substitution  K  -  K’  in  the  integrand. 

For  every  value  of  K.,  the  integrand  of  Eq .  (3.10)  represents  a  plane 
wave  with  wave  vector  K,  and  because  of  the  identity 

a  A  * 

K*(k  £(w)  -  (m(w).K)K)  -  0  ,  (3.11) 

every  wave  is  transverse.  The  z -component  of  the  wave  vector  K  equals  ^  from 
Eq.  (3.3),  which  is  either  positive  or  positive- imaginary,  corresponding, 
respectively,  to  a  wave  which  travels  in  the  positive  z-direction  or  to  a  wave 
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which  decays  exponentially  in  amplitude  in  the  same  direction  (evanescent 
wave).  For  the  field  in  the  region  0  <  2  <  h,  the  wave  vector  is  K'  for  a 
given  K|| ,  and  the  wave  is  again  either  travelling  or  evanescent.  The  various 
occuring  waves  are  illustrated  pictorially  in  Fig.  1. 


These  unit  vectors  are  normalized  according  to  e.»e^  “  1  f°r  any  subscript  i, 
and  they  are  perpendicular  to  their  corresponding  wave  vector.  Notice  that 
the  p-polarization  vectors  can  have  an  imaginary  z-component.  The  directions 

A 

of  the  various  unit  vectors  are  depicted  in  Fig.  1.  The  field  Ep(r,w)  In  z  > 
h  can  then  be  decomposed  as 
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(r,y)  -  - / 

Sff“£  c“  - 
o  a 


-P  ~ 


,2,,  1  iK*  (r-h)  . A  ,  ,  . 

d  Km  ~  e  -  -  -  (M(w)*e.  )e 

—  '  "y  —  ~J\.0  “'!\.£7 


(4.4) 


vi 


th  a  -  s,p,  and  with  K  -  K'  ve  obtain  the  field  for  0  <  z  <  h. 


V.  REFLECTED  FIELDS 

In  order  to  simplify  the  notation,  we  introduce  the  two  quantum 
operators 


.  2  -my 

iv  e  ,  .  . 

I  'mv)*e.  ) 

o  -  -  -  ~^a 

3*  £  C  1 

o 


(5.1) 


.  2  ih7 

iv  e  ,  ,  .  . 

“K'<r  ”  2  2  ‘-K'cr 

£  C  7 
o 


(5.2) 


Then  the  particular  solution  can  be  written  as 

E  (r , w)  -  )  [  d2K„  E,,  e.,  e1-*-  .  for  2  >  h  ,  (5.3) 

-p  -  L  J  -:!  Kt7-Kt7 

<7 

and  with  K  -  K'  we  find  the  field  in  the  region  0  <  z  <  h.  Expression  (5.3) 
with  K  -*  K’  is  a  superposition  of  polarized  plane  waves,  where  serves  as 

the  amplitude,  including  the  proper  phase.  Since  Maxwell's  equations  are 
linear,  we  can  calculate  the  reflected  field  on  a  per-wave  basis.  In  the 
region  0  <  z  <  h,  each  plane  wave  has  the  form  Zv,  e„,  exp(iK'»r),  and  these 
waves  are  the  incident  waves  on  the  medium  (Fig.  1),  which  have  to  be  matched 
across  the  boundary  z  -  0  to  a  solution  in  z  <  0. 
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.v'hen  a  plane  wave  with  wave  vector  K'  ,  frequency  u  >  0  and  polarization 

a  is  incident  upon  the  medium,  then  the  reflected  radiation  consists  again  of 

a-polarized  plane  waves,  but  with  different  wave  vectors,  in  general.  Due  to 

the  four-wave  mixing,  the  frequency  of  the  reflected  wave  can  also  be 

different.  When  the  two  pumps  have  frequency  u>  >  0,  and  the  incident  wave  has 
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frequency  w  >  0,  then  it  can  be  shown  that  the  FWM  process  couples  the  field 
with  frequency  tw  to  the  field  with  frequency 

"  «*»  *  ,  (5.4) 

which  is  negative.  Therefore,  reflected  radiation  can  be  generated  in  the 

medium,  both  at  frequency  u  and  at  frequency  w' ,  but  not  at  other  frequencies. 

At  the  interface  z  -  0,  the  spatial  part  of  the  incident  wave  reduces  to 

exp(iKi|»r).  It  is  only  possible  to  match  this  wave  to  other  waves  for  all  r 

simultaneously  when  the  other  waves  have  the  same  spatial  dependence. 

Therefore,  all  waves  inside  and  outside  the  medium  must  have  a  wave  vector 

with  the  same  parallel  component  Ki|.  For  the  wave  vector  of  the  reflected 

wave  at  rrequency  w,  we  can  write  K'  -  Ku  +  K’  e  ,  for  which  it  must  hold 

-r  -;j  r,z  -z 

2  2  2  2  2 

that  (K^.  z)  -u>/c  -  Ki|  -  7  .  The  solution  ^  -  -7  corresponds  to  the 

incident  wave,  and  the  only  other  possibility  is  -  +7 .  But  this  is  the 

z-component  of  the  complementary  wave  K,  so  that  we  find  -  K.  Figure  1 
shows  that  this  wave  is  just  the  ordinary  specular  wave,  alt  hough  its 
generation  is  affected  by  the  nonlinear  interaction.  For  the  wave  at 
frequency  w'  we  wr  te 


K’  “  Km  +  K'  e 
-pc  pc,z  -z 


(5.5) 
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where 


-v(^'/cr  -  Km 

.1 

(5.6) 

i7K^  -  (u'/c)*" 


In  general,  ve  shall  have  lea'  I  =»  leal  ,  which  gives  K'  =  -y  for  Che  case  of  a 

1  °  pc.z 

cravelling  wave,  and  thereby  K'  =  K' .  This  shows  Chat  che  incident  wave  and 

-  -pc 

the  pc-wave  are  approximately  counterpropagating ,  as  is  required  for  a  time- 

-  ~  V  0  T  S  £  G  V  3  V  6  . 

'.'ext  we  have  to  choose  a  phase  convention  for  the  unit  polarization 
vectors  of  the  r-wave  and  the  pc-wave.  Since  K^.  -  K,  we  can  simply  take 


e., 

-x 


a  -  s,p 


(5.7) 


tor  the  r-wave,  where  the  right-hand  side  is  given  bv  Eqs.  (4.1)  and  (4.2).  A 
convenient  choice  for  the  pc-wave  is 


s 


x  e 
~z 


(5.8) 


e„ ,  -  -  ~  K*  x  e., 

~fv  p  wJ  -pc  -ts’  s 

-pcr  r  -DC 


(5.9) 


which  are  both  normalized  according  to  e.«e.  -  1,  and  are  both  perpendicular 

to  K'  . 

-pc 
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For  a  given  incident  field  of  the  form  (5.3)  with  K  -»  K' ,  the  reflected 
field  (homogeneous  solution)  has  two  components.  The  r-wave  at  frequency  u> 
attains  the  form 


E 

-r 


(r,w) 


(5.10) 


and  for  the  pc-wave  we  can  write 


E 

-pc 


(r,«') 


EK'aFK'c2K'  a6 
-pc 


iK'  •  r 
-pc  - 


(5.11) 


Here,  R„ ,  and  P  are  the  Fresnel  coefficients  for  reflection  of  the  r-wave 
rC  a  rC  a 

and  the  pc-wave,  respectively.  They  depend  in  a  complicated  way  on  the 
properties  of  the  medium  and  the  details  of  the  FWM,  like  the  dielectric 
constant,  the  third-order  susceptibility,  the  pump  polarization,  etc. 


VI .  TIME  DOMAIN 

Of  practical  interest  is  the  electric  field  as  a  function  of  time.  It 
is  sufficient  to  evaluate  the  positive-  or  negative- frequency  part  only  since 
these  components  determine  the  field  completely.  For  the  particular  solution 
Ep(r,w)  we  have  Eq .  (5.3)  for  z  >  h,  which  holds  for  u  >  0  only.  Therefore, 
the  posi tive - frequency  part  of  E  (r,t)  follows  from  Eq  .  (2.6),  and  is 

-p  - 

explicitly 


V£-c> 


c+)  l  r,  .i«t \  r  ,2„  _ 

"  2 ,]0due  2J  d-il  V~K* 


iKT 
e  -  - 


(6.1) 
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for  z  >  h,  and  with  K  -  K'  we  obtain  the  field  in  0  <  z  <  h.  In  a  similar 
vav,  rhe  superposition  of  all  specular  waves  has  the  representation 


Er(r,t) 


(+) 


fell 


dw  e 


-  ltdt 


z  I  “2-;i 


Ka 


iK*r 
e  -  - 


(6.2) 


For  the  phase -conjugated  wave  we  have  expression  (5.11)  in  which  the  frequency 
1  is  negative.  Consequently,  a  Fourier  synthesis  yields  the  negative- 
frequency  part  of  the  field.  We  obtain 


(ri  -  ) 

-pc  - 


^.T 


du>'  e 


-  ij  t 


U 


o 

d^Kii 


Z„,  ?,.,  e.,,  e 
K  a  K  a-Pv  a 
-pc 


: 

-pc 


(6.3) 


VII.  ASYHPTQTIC  EXPANSION 

Although  expressions  (6.1) -(6. 3)  give  the  total  source  field  in  z  >  0  as 
a  function  of  time,  the  appearance  of  the  integrals  over  the  frequency  and  the 
parallel  components  of  the  wave  vectors  is  cumbersome.  Fortunately,  for  the 
studv  of  fluorescent  emission  we  do  not  need  the  exact  solution  for  the 
radiation  field  at  every  point  r  in  the  half-space  z  >  0.  Only  the  value  of 
E  (r,t)  for  r  -  |r|  -•  «  is  of  relevance.  In  this  section  we  evaluate  the 
asymptotic  behavior  for  r  -•  ® ,  and  in  the  next  two  sections  we  work  out  the 
formal  results  in  order  to  obtain  more  managable  and  transparent  expressions. 

Suppose  a  detector  is  located  at  position  r.  This  point  will  be 

represented  by  its  spherical  coordinates  (r with  respect  to  the  z-axis, 

and  we  are  interested  in  the  value  of  the  radiation  field  for  r  -»  «  with  (6 ,4>) 

fixed,  and  for  z  >  0  only.  The  integrals  in  Eqs.  (6.1)-(6.3)  have  the  form 

I  d2K„  exp ( i*»  r )  g(K,i ) ,  with  k  -  K  or  *  -  K'  .  With  the  method  of  stationary 
J  ~  il  -  -  ~:l  -  -  -pc 
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phase,  we  can  derive  an  asymptotic  expansion  for  integrals  of  this  kind. 
For  the  p-wave  and  r-wave  we  obtain 


:0(£-C) 


icost?  r 

■~Jo 


11  ''‘"<M/C)  “  ) 

a  ”  ' 


,  (7.1) 


E  (r.t) 
~r  - 


(7.2) 


—  -  1  +  2-~- 


(7.3) 


^he  asymptotic  expansion  of  the  pc-wave  becomes 


-pc  -  rc 


l“'<c-r/c>  “-JiVA-A-  A,-*,.,, 

~  -  -pc  ~r 


(7.4) 


The  parameter  p  accounts  for  the  mismatch  between  u>  (or  -u>’)  and  the 
frequency  w. 


VIII. 


In  the  formal  expansions  (7.1),  (7.2)  and  (7.4)  we  still  have  to 
evaluate  the  factors  in  square  brackets  for  the  indicated  values  of  K|  .  In 
this  section  we  consider  the  (quantum)  amplitude  factors  and  Ej^,  as 
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by  Eqs .  (5,1)  and  (5.2),  respectively.  These  factors  contain  the  parameter  7. 

O  O  '' 

which  depends  on  -  a~+  3"  according  to  Eq.  (3.3).  Ve  find 


2  2 

k  sin  8  ,  for  the  p-  and  r-wave 

(8.1) 

2  2  2 

p  k  sin  8  ,  for  the  pc -wave 


which  gives 


f  kcos»  for  the  p-  and  r-wave 


[  kyi  -  p“s in~8  ,  for  the  pc-wave 


(8.2) 


Substitution  of  E  and  E  into  the  expansions  (7.1),  (7.2)  and  (7.4)  then 
gives  for  the  three  fields 


:  (r,t)(+)  -  — ■  ■*  ~  j* 

p  rc  J0 

o 


,  -ia>(t-r/c+r)  2 

e  cj 


*  /  :-Ka(^(“)‘SKo)]K|,-k(er)||  ' 


:r(£.t)(  +  )  -  — 2  1  ■  J" 

8 we  rc  J0 


.  -iw(t-r/c-r)  2 
du  e  w 


\  [RK'ff2Kff(tf(w),SK'<7;)lK.-k(er)| 


(8.3) 


X 


(8.4) 
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E 

-pc 


(r.t) 


(-) 


1 


2  2 
Sjt  c  rc 
o 


dw'e'i^' (c-r/c+r^/p) 


UX j’ 


*/  ;PK'a-K'  £7(^(w)*-K'a)  ]K,|  — pk(e  ), 
g  ~  ~PC  ~  “ll  -r 


(8.5) 


where  we  introduced  the  parameters 


/I  -  o  sin~j 
cost) 


T  -  —  cost) 
c 


(8.6) 


(8.7) 


At  this  stage  it  is  convenient,  although  not  necessary,  to  make  a  slight 

approximation  in  the  expression  for  the  phase  -  conj  ugated  wave.  For  a  given 

frequency  w' ,  the  integrand  in  Eq .  (8.5)  is  proportional  to  the  Fresnel 
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reflection  coefficient  P^,^,  which  depends  on  w' .  It  is  well  known  that 

P  ^  assumes  only  a  finite  value  when  the  frequency  of  the  incident  wave  is 

sufficiently  close  to  the  pump  frequency  w.  Since  the  integrand  in  Eq .  (8.5) 

is  proportional  to  P  ,  we  can  effectively  set  p  =  1  in  any  factor  that 

multiplies  P^,  .  Similarly,  we  can  set  £  -  1  and  w'  -  -u>  in  the  integrand, 

and  also  e„,  -  e„,  as  can  be  verified  by  inspection  of  the  definition  of 

-K  a  -K  <7 
-pc 

these  unit  vectors.  Then  Eq.  (8.5)  reduces  to 


ECr.t) 
-pc  - 


(-) 


,  ,  -iu>'  (t-r/c+r)  2 

d4*>'  e  w 


c  rc  J  -« 
o 


)  ^PK’a-K’ff(£(w)'-K'c7)  ^K|— k(er). 


x 


(8.8) 
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which  is  a  great  simplification.  Notice  that  we  have  set  p  -  1  only  in  the 
overall  factors,  but  the  full  resonant  frequency  dependence  is  retained, 

as  it  is  incorporated  in  P„,  . 

K  o 


IX.  POLARIZATION-  VECTORS ■  FRESNEL  COEFFICIENTS  AND  THE  MIRROR  DIPOLE 

Next  we  have  to  evaluate  the  polarization  vectors  e  and  e  ^ ^ ,  as  they 

are  defined  by  Eqs.  (4.1)-(4.3),  and  for  the  indicated  values  of  K >|  in  Eqs . 

8.3),  (8.4)  and  (8.8) .  The  parameter  7  in  Eqs.  (4.1)*(4.3)  equals  7  —  kcosfl 

m  ail  cases,  as  follows  from  Eq .  (8.2)  with  p  ■*  1 .  We  wish  to  express  the 

polarization  vectors  in  the  standard  spherical  unit  vectors  e.  and  e,  for  a 

-?  — P 

given  observation  angle  (S  .3)  . 

For  the  particular  solution  (p-wave)  from  Eq .  (8.3)  we  need  e,.  for  K«  - 

— K  o  —j| 

^(®r)  !  1  both  for  a  -  s  and  a  —  p.  We  find 


'-Ks  Kn-k(et_)  ,1  -4> 


(9.1) 


'-KpiK|!-k(er)  |  "  -? 


(9.2) 


as  could  be  expected  from  a  comparison  between  the  phase  conventions  for  ev 

-NS 

and  e,  (Fig.  1),  and  the  definitions  of  e,  and  e.,  respectively.  The  field 

~<p  -o  r 

then  attains  the  form 


E  (r ,  t) (  +  )  -  — -  f° 

-p-  7  2  0 

8tt  c  re4" 


,  -  iu>(  t- r/c+r )  2 

at*?  e  w 


X  { (g(w)*etf)etf  +  (£<w)«e^)e  )  , 


(9.3) 


however ,  that 


fluorescence,  which  is  of  course  unphysical.  Ic  can  be  shown, ^ 

when  we  allow  for  a  finite  but  very  small  detuning  -  u  j*  0 ,  the  Fresnel 

coefficients  remain  finite,  with  a  magnitude  on  the  order  of  unity.  This 

renders  finite  values  for  the  parameters  c  ,|  and  cj_,  and  thereby  for  the 

emission  rate.  For  an  extremely  close  resonance  between  u>  and  w  we  can  no 

o 

longer  replace  P^(<j,cos0)  in  Eq.  (9.15)  by  ,cos$)  ,  and  we  have  to  take 

into  account  the  frequency  dependence  of  P  (u,cos 6)  over  the  width  of  the 

a 

atomic  emission  line. 

»his  situation  is  analogous  to  the  problem  of  resonance  fluorescence  by 

a  two-^evel  atom  in  a  monochromatic  laser  field.  There,  the  spectral 

distribution  of  the  fluorescence  consists  of  a  6-function  at  the  laser 

frequency,  superposed  on  a  smooth  background,  and  in  the  limit  where  u>  equals 

o 

the  laser  frequency  the  emission  rate  is  infinite.  In  any  case,  the  example 
of  this  section  shows  that  the  situation  of  perfect  resonance  has  to  be 
considered  with  caution. 


XV.  SCHRODINGER  PTCTirRF. 

The  time  dependence  of  the  emission  rate  dW/dt  in  Eq .  (13.6)  is  governed 
by  the  time  evolution  of  the  Heisenberg  operator  .  A  more  transparent 
representation  of  dW/dt  can  be  obtained  by  a  transformation  to  the  Schrddinger 
picture.  When  we  take  the  Schrddinger  and  Heisenberg  picture  to  coincide  at 
t  -  0 ,  then  we  have 

ff(°)  -  if  -  M(  +  )  +  M(  *  ^  .  (15.1) 


with 


20 


J 


and  notice  chat  the  factor  in  brackets  is  just  the  transverse  component  of 
^w)  with  respect  to  che  propagation  direction 

For  the  r-wave  from  Eq.  (8.4)  we  need  again  e^  at  Kij  -  k(er)|  for  the 
unit  polarization  vectors  which  correspond  to  a  propagating  wave  with  wave 
vector  K  in  the  observation  direction  In  addition,  we  need  at  the 
same  value  of  Kij  ,  and  these  unit  vectors  account  for  the  polarization  of  an 
incident  wave  on  the  surface  z  -  0,  which  is  subsequently  emitted  in  the 
specular  direction  (Fig.  1).  We  find 


t 


(9.4) 


•  —  K '  p '  K.i  — k( £r ) 


-e,  -2sini?e 
-9  -z 


(9.5) 


where  the  term  -Isintfe^  represents  the  change  in  direction  between  K'  and  K  at 
the  moment  of  reflection  at  the  interface  z  -  0.  For  s -polarized  waves  there 
is  r.o  such  correction,  since  both  the  incident  wave  and  the  specular  have  the 
same  polarization  vector  ( Eq .  (4.1)).  The  specular  field  becomes 


Er(r,t) 


(  +  ) 


1 

q  2  2 

Sit  (  rc 
o 


.  -iw(t-r/c-r) 
ou>  e 


2 

w 


x 


lRK'p(-r',>)'  ( '2sin&tz^tg 


(9.6) 


where  the  Fresnel  coefficients  R^,,^  pertain  to  che  ’flection  coefficients  of 
incident  waves  with  wave  vectors  K'  ,  which  are  specularly  reflected  in  the 
observation  direction  (6,#).  The  coefficients  depend  on  the  frequency  w 


I 


and  angle  of  incidence  6^.  Since  I K  equals  the  observation  angle  8,  we  might 

as  well  label  R.,,  with  8,  rather  than  6..  Therefore,  we  can  write 
K  a  t 


Rjj,^  ”  Rct(w,cos 8)  ,  a  -  s.p 


(9.7) 


A  convenient  concept  for  specular  reflection  is  the  mirror  dipole.  If  we 
decompose  m  (or  m)  as 


a  -  mi  +  M.| 


(9.8) 


where  the  subscripts  1  and  '|  refer  to  the  surface 
dipole  is  defined  as 


z  -  0,  then  the  mirror 


M  **  Mj.  - 


(9.9) 


Combining  everything  then  yields  for  the  specular  field 


E  (r,t)(  +  )  -  -2-L7  f 
L  a-4-*  J n 


- iu(t-r/c-r )  2 

du  e  w 


8 jt*" e  rc  J0 
o 


x  (R  (<4>,cos9)  (p1  (w)»e^)e# 


Rs(w,cos0)(m' (w)*e)e^) 


(9.10) 


The  minus  sign  in  front  of  R  is  a  consequence  of  the  phase  convention  of  the 


polarization  vectors  for  s -waves. 


For  the  phase -conjugated  field  we  have  to  calculate  e^,ff  at  K« 
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-k(e  )  .  Both  unit  vectors  in  Eq.  (8.8)  appear  with  wave  vector  K*  ,  which 
reflects  the  fact  that  both  the  incident  plane  wave  and  the  pc -wave  have  the 
same  wave  vector  K'  .  The  additional  minus  sign  in  K.j ,  as  compared  to  Kg  for 
the  p-wave  and  the  r-wave,  signifies  that  a  phase-conjugated  plane  wave 
travels  in  the  direction  opposite  to  its  wave  vector.  The  polarization 
vectors  are 


•-sC's^K..  — k(er),|  " 


(9.11) 


'  -V. ' p ‘ K  .--k(e  i  . 
~  'I  -r'  I 


=  e  . 


(9.12) 


snd  the  field  is  found  to  be 


-pc  - 


(*) 


r 

,  .2  J- 


.  ,  -iw' (t-r/c+r)  2 

dw'  e  '  u 


8ir  (  rc  J  -« 
o 


X  ( Pp(w,cos^)(/<(w)*e^)e^ 


+  ?s (w, cos^ ) (^(«) »e^)e^ ) 


(9.13) 


•;here  we  introduced  the  notation 


P  (w.costf)  -  Pv,  ,  am  s , p 

Q  r*  O 


(9.14) 
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A  complication  with  expression  (9.13)  is  that  the  integration  runs  over 
whereas  the  factor  in  brackets  depends  on  w.  With  the  relation  w'  -  w  - 
2w  we  can  write  the  field  in  the  alternative  form 


r 

dpc 


(r,t) 


(-) 


1  2iw(t-r/c+r)  P  ,  - iw( t- r/c+r )  2 

8*2e  rc2  6  JO  * 

o 


x  ( P^(w, costf) (m(w) *e^ )e^ 

A 

+  P  (w,cosi/)(/j(w)«e,)e,)  ,  (9.15) 

S  -  -.5 


•..here  we  have  used  that  P^(w.costf)  is  only  nonzero  in  a  small  frequency  band 
around  w . 


X.  NARROW- BAND  EXCITATION 


The  Fresnel  coefficients  R  (w.costf)  and  P  (w.costf)  depend  in  general  in 

a  a 

a  complicated  way  on  the  frequency  w.  For  R  this  is  mainly  brought  about  by 


the  variation  of  the  dielectric  constant  with  w,  but  P  has  an  additional 


a 

geometric  frequency  dependence.  The  value  of  P^  has  a  sharp  peak  around  w  =« 
w,  and  the  relative  frequency  width  of  the  response  is  on  the  order  of  |/’-l|  - 
7.  The  absolute  frequency  width  then  becomes  |w-ui|  -  yw,  and  even  for  very 
small  interaction  parameters  7  this  can  still  be  very  large  due  to  the 
multiplication  by  w.  We  now  assume  that  the  dipole  radiation  is  nearly 
monochromatic,  compared  to  the  frequency  width  of  the  PC.  With  <jq  as  the 
central  frequency  of  the  dipole  field  and  with  Aw  as  its  typical  width,  we 
impose  the  restriction  Aw  «  yw  on  the  exciting  field.  However,  this  does  not 
imply  that  wq  has  to  be  in  close  resonance  with  w,  but  only  that  the  spectral 
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width  of  the  dipole  radiation  is  small  compared  to  70.  For  an  atomic 
transition  between  (nearly)  degenerate  states  which  are  separated  by  wq,  this 
condition  is  very  easily  met. 

Under  the  above  condition  we  can  replace  R  (w.cosS)  and  P  (w.cosfl)  by 

a  o 

2 

their  values  at  u>  -  w  .  Furthermore,  the  factor  w  in  the  integrands  can  be 

o 

2 

taken  outside  the  integrals  as  The  only  remaining  frequency  dependence  of 

the  integrands  enters  as  m(u)  or  its  mirror  image,  and  the  w- integrals  can  be 
performed  easily.  We  find  for  the  three  fields 
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X  *  Pp^0.costf)(if(t.r/c+r)(+).e<?)e5 

+  P  (w  , cos0 ) (ji( t-r/c+r ) ^ *e , )e . )  ,  (10.3) 

so  —  *•  <p  ~-<p 

in  terms  of  the  positive- frequency  part  n(t)^  of  the  dipole  moment  operator. 
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XI.  RETARDATION 

The  three  fields,  evaluated  at  the  position  r  of  the  detector  and  at 
time  t,  are  determined  by  the  instantaneous  value  of  the  dipole  moment  at  an 
earlier  time  t  -  r/c  ±  t.  Since  photons  travel  with  the  speed  of  light  c,  the 
delay  time  r/c  equals  the  propagation  time  of  a  photon  when  it  would  travel 
from  the  origin  of  the  coordinate  system  to  the  detector.  For  both  the  p-wave 
and  the  pc -wave,  this  retardation  time  is  reduced  by  an  amount  r  —  (h/c)costf 
for  a  given  observation  angle  S,  and  Fig.  2  illustrates  that  this  implies  that 
the  photons  are  emitted  from  the  site  of  the  dipole,  and  directly  into  the 
direction  of  the  detector.  The  p-wave  corresponds  to  emitted  radiation  by  the 
dipole  into  the  region  z  >  h  without  interference  from  the  medium  (Sec.  IV), 
which  makes  this  picture  consistent.  The  phase -conjugated  wave,  however,  is  a 
reflected  field  by  the  surface,  but  nevertheless  the  photons  seem  to  emanate 
from  the  site  of  the  dipole  and  travel  directly  into  the  direction  of  the 
detector,  in  view  of  their  retardation  time.  The  specular  r-wave  has  a 
retardation  of  2r,  as  compared  to  the  p-wave  and  the  pc-wave,  and  it  is  shown 
in  Fig.  3  that  this  time  delay  accounts  for  the  difference  in  travel  time 
between  a  directly-emitted  photon  and  a  photon  which  is  first  reflected  by  the 
surface  z  -  0.  The  figure  also  illustrates  that  upon  reflection  at  point  A 
the  path  of  a  photon  obeys  Snell's  law,  and  that  a  photon  appears  to  be 
emitted  by  a  mirror  dipole  n'  which  is  a  distance  h  below  the  surface.  We 
conclude  that  the  emission  of  both  p-  and  r-photons  is  in  strong  analogy  with 
classical  ray  optics,  and  that  the  quantum  nature  of  the  radiation  affects 
neither  the  geometry  nor  the  interpretation  of  the  light  emission.  This  is  in 
contrast  to  the  emission  of  the  pc-wave.  If  this  wave  would  be  a  result  of 
the  emission  of  a  photon  by  the  dipole  into  the  direction  of  the  surface,  and 
a  subsequent  reflection  as  a  phase -conjugated  image,  then  the  delay  time  would 
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be  2h/cos0,  as  compared  to  direct  emission.  Since  there  is  no  such 
retardation,  we  conclude  that  the  emission  of  phase-conjugated  fluorescence  is 
more  complicated.  We  shall  show  in  Sec.  XVII  how  the  mechanism  of  pc-wave 
emission  can  be  understood. 

XII.  TOTAL  SOURCE  FIELD 

In  this  section  we  construct  the  total  source  field  E  ,  as  it  is 

-  s 

measured  by  the  detector.  We  shall  suppress  the  overall  time  delay  with  r/c. 
Then,  for  an  atom  which  is  not  more  than  a  few  wavelengths  away  from  the 
surface  z  -  0,  the  time  delay  r  is  only  a  few  optical  cycles,  and  we  can 
safely  replace  the  time  evolution  of  /i(t)  by  its  free  evolution  over  this 
small  time  interval.  This  amounts  to  the  approximation 

-  YU)  7 

u(t =<  e  0  a(t)^  ,  (12.1) 

which  is  perfectly  justified  at  optical  frequencies.  We  then  obtain  for  the 
fields 
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in  which  all  dipole  moments  appear  with  the  same  time  argument  t. 

Subsequently,  we  take  the  Hermitian  conjugate  of  Eq.  (12.4)  which  gives 
us  E  ■'.r,t)^+^  according  to  Eq .  (2.7),  and  we  use  (/j(t)^  +  ^l^  -  /i(t)^  ^  ■  Then 
we  add  the  three  fields,  which  gives  us  for  the  positive -frequency  part  of  the 
source  field 
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•here  we  introduced  the  two  operators 
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Without  the  presence  of  the  medium  both  operators  reduce  to  M.  -  M,  -  t)^+\ 

—  P  ~<P  — 

and  conversely  we  can  interpret  M.  and  M,  as  two  components  of  the  positive- 

-/)  ~<p 

frequency  part  of  an  effective  dipole  operator  which  takes  into  account  the 
interaction  with  the  medium.  Notice  that  both  operators  contain  only  positive 
frequencies,  due  to  the  appearance  of  the  factors  exp(-2iut). 


XIII.  INTENSITY 

When  a  detector  at  a  distance  r  and  under  solid  angle  0  -  (6  ,<j>)  measures 

an  electric  field  Eg,  then  the  detected  energy  at  time  t  per  unit  of  time  and 

39 

per  unit  solid  angle  is  in  general  given  by 


3  tan 


W  -  2eQcr2  <Eg(r,t)('')*Es(r,t)^> 


(13.1) 


provided  that  the  efficiency  is  100%.  Here  the  angle  brackets  indicate  a 
quantum  expectation  value.  The  emitted  power  into  the  half-space  z  >  0  is 
defined  by 


dW 

dt 
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With  expression  (12.5)  for  the  source  field,  we  obtain  for  the  intensity 
distribution 
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+  (M , ( t , costf ) 1 «e , ) (M , ( t , costf ) *e  )> 
-< P  ~<P  -v 


(13.3) 


and  substitution  of  the  explicit  forms  of  and  then  gives  the  angular 
radiation  pattern. 

Both  operators  and  M ,  depend  only  on  the  angle  9  ,  and  therefore  the 

—  V  -~<p 

2 

^-dependence  in  3  U/3tdn  is  purely  geometrical,  which  reflects  the  invariance 
of  the  system  for  rotations  about  the  z-axis.  Integration  over  4>  and  6  then 
yields  for  the  emitted  power 
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where  the  overall  parameter  is  given  by 
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and  the  parameter  functions  are 
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(13.9) 


Here,  u  signifies  the  cosine  of  the  angle  of  incidence  of  a  plane  wave,  and 
the  integrals  represent  the  superposition  of  plane  waves.  The  parameter 
functions  bi|,  b^,  cu  and  c_j_  are  real  and  positive,  and  they  are  determined  by 
the  Fresnel  reflection  coefficients.  These  four  parameters  incorporate 
entirely  all  the  properties  of  the  medium,  like  its  dielectric  constant, 
nonlinear  interaction  parameter  and  geometry  (for  instance,  the  layer 
thickness),  and  they  are  independent  of  the  dipole  moment  /j.  Therefore,  the 
representation  (13.4)  for  the  emission  rate  makes  a  clear  separation  between 
material  properties  and  the  dynamic  evolution  of  the  atom,  which  is 
incorporated  in  the  time  dependence  of  the  dipole  moment  n(t) .  This  result  is 
reminiscent  of  the  general  form  of  the  spontaneous -decay  operator  for  an  atom 
near  a  linear  medium,  which  can  also  be  expressed  in  terms  of  the  parameter 
functions  bi|  and  bj_  and  dipo le -moment  expectation  values.  The  connection 

ii 

is,  of  course,  that  spontaneous  emission  and  spontaneous  decay  are  intimately 
related. 

The  parameter  $  in  Eqs .  (13.6)  and  (13.7)  is  defined  as 


&  -  2uQh/c  ,  (13.10) 

which  equals  4?r  times  the  distance  h,  measured  in  optical  wavelengths.  This 
parameter  enters  through  the  cos $• dependence  of  the  delay  time  r,  and  the 
factor  exp(i/3u)  accounts  for  the  proper  phase  relation  between  the  interfering 
p-waves  and  r-waves.  There  is  no  such  factor  in  the  expressions  for  C|  and  c_l 
(which  describe  the  pc-wave  emission),  and  this  indicates  again  that  the 
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p  h  as  e  -  c  o  n  j  u  ga  c  e  d  fluorescence  is  of  a  different  nature  than  ordinary 
fluorescence  (p-waves  and  r-waves).  Especially,  the  pc-waves  do  not  interfere 
with  the  other  components  of  the  emitted  radiation.  Also  notice  that  the 
parameters  bj|  and  bj_  depend  on  the  normal  distance  between  the  atom  and  the 
surface  (through  0) ,  but  that  Ci|  and  c^  are  independent  of  h.  It  can  be  shown 
that  this  is  a  consequence  of  approximation  (12.1),  which  limits  the  range  of 
h  to  a  few  optical  wavelengths.  In  physical  terms,  the  atom  does  not  decay  on 
a  time  scale  r. 

XIV.  PARAMETER  FUNCTIONS 

In  order  to  shed  some  light  on  the  significance  of  the  parameter 
functions,  we  consider  some  examples. 

A.  Transparent  medium 

When  the  dielectric  constant  e  of  the  medium  equals  unity  and  the 
nonlinear  interaction  parameter  7  equals  zero,  then  there  is  effectively  no 
medium  at  all,  which  gives 


R  -  R  -  P  -  P  -0  (14.1) 

s  p  s  p 


for  the  Fresnel  coefficients.  Then  the  parameter  functions  are  found  to  be 


Consequently,  a  deviation  of  b||  or  bj_  from  unity  reflects  the  presence  of 
specular  waves,  and  a  nonvanishing  c  ^  or  cj_  indicates  a  phase  * conj ugate 
signal.  The  emitted  power  becomes 


32 


Jjf  -  a  <M(t)(').y(c)^  +  )>  .  (14.3) 

B .  Perfect  conductor 

For  the  case  of  a  perfect  -  conducting  medium  (mirror)  the  parameter 
functions  can  readily  be  evaluated.  The  Fresnel  coefficients  are 


P  -  P  -  0 
s  p 


(14.4) 


which  gives  of  course 


ci|  -  cj_  -  0 


For  the  coefficients  of  linear  reflection  we  find 


b.  -  2  -  3(^  +  .  sinfi 

11  0  p2  p2 


bx  -  2  -  6(“^  -  ainflj 

r  ,3 


which  is  a  well-known  result. 


(14.5) 


(14.6) 


(14.7) 


C.  Ideal  PC 

For  a  plane  wave  at  normal  incidence,  the  distinction  between  s-waves 

and  p-waves  disappears,  and  we  have  P  -  P  .  When  most  of  the  incident 

s  p 

radiation  is  normal  to  the  surface,  we  can  make  the  approximation  Pg  “  Pp  and 
take  the  value  u  -  1.  For  an  ideal  PC  the  reflection  coefficient  is 


A 
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independent  of  the  polarization  and  the  angle  of  incidence,  and  there  is  no 
specular  reflection.  Then  the  Fresnel  coefficients  are 


R-R-0  ,  P-P-P  , 

s  p  s  p 


(14.8) 


and  in  analogy  with  the  values  for  a  perfect  conductor,  Eq.  (14.4),  such  a 
device  is  sometimes  called  a  phase-conjugating  mirror.  We  find 


b,.  -  bx  -  1  ,  c —  cj_  —  )  P  | 


(14.9) 


and  the  emission  rate  becomes 


f;  -  a(<M(C)(')-/f(t)(  +  )>  +  1  P  j  2  <tt(  t)(  +  )-/i(t)(‘)>) 


(14.10) 


D .  Transparent  PC  on  resonance 

A  transparent  medium  is  defined  as  having  a  unit  dielectric  constant, 
which  gives  for  the  parameters  of  the  directly-emitted  and  specular  waves 


R  —  R  -  0  ,  bii-bi-1 

s  p  ,|  x 


(14.11) 


If  we  assume  that  the  atomic  transition  frequency  is  in  close  resonance 
with  the  pump  fields,  e.g., 


|  u>Q  -  w  |  «  7w  , 


(14.12) 


then  the  Fresnel  reflection  coefficients  take  the  form 


iPs(«o,u)|2  -  tan2(J)  , 
|Pp(«o,u)  |2  -  tan2(JJ(3-2u2)) 
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(14.13) 


(14.14) 


in  terms  of  the  dimensionless  parameter 


(14.15) 


where  A  is  the  layer  thickness  of  the  medium.  For  almost  perpendicular 

2  o  2 

ir.cider.ce  these  extress  ions  reduce  to  !  P  I  »  I  ?  !  “  =■  tan'r 7.  and  this  behavior 

s  '  ?' 
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has  been  confirmed  qualitatively  bv  experiment.  ’  It  as  also  well  known 

2  2 

that  for  certain  values  of  t)  (or  A  and  u  )  the  magnitude  of  |P  |  or  |P  |  can 

o  s  ^ 

become  infinite,  as  follows  from  Eqs.  (14.13)  and  (14.14).  This  phenomenon  is 
termed  se  1  f  -  osc  i  1 1  a  t  ion  and  has  been  observed  in  experiment  also.  For  the 
present  problem,  however,  the  cosine  of  the  angle  of  incidence  is  a  variable, 
rather  than  the  parameter  rj ,  and  it  is  seen  from  Eq.  (14.13)  that  for 


cos*  -  — ^  (14.16) 

(,n-t-’5)7r 

with  n  integer,  the  value  of  ! ?  I  is  also  infinite.  In  fact,  condition 

's' 

(14.16)  predicts  an  infinite  series  of  values  of  $  in  the  range  (0“,90°)  for 

2  2 
which  J  P  s  |  diverges.  In  a  similar  way,  the  reflectivity  |Pp|  has  a  series 

of  resonances  for  certain  angles  of  incidence. 

When  we  substitute  the  expressions  (14.13)  and  (14.14)  into  Eqs.  (13.8) 

and  (13.9),  then  it  is  easy  to  verify  that  the  integrals  over  u  for  c <|  and  cj_ 

diverge.  This  yields  an  infinite  emission  rate  for  the  phase-conjugated 
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fluorescence,  which  is  of  course  unphysical.  It  can  be  shown,  however,  thac 
when  we  allow  for  a  f inice  but  very  small  detuning  -  w  *  0,  che  Fresnel 
coefficients  remain  finite,  with  a  magnitude  on  the  order  of  unity.  This 

renders  finite  values  for  the  parameters  c,|  and  cj_ ,  and  thereby  for  the 

ii 

mission  rate.  For  an  extremely  close  resonance  between  and  w  we  can  no 
longer  replace  P^(u,co s$)  in  Eq.  (9.15)  by  P^(wq,cos0),  and  we  have  to  take 
into  account  the  frequency  dependence  of  P^(w.costf)  over  the  width  of  the 
atomic  emission  line. 

This  situation  is  analogous  to  the  problem  of  resonance  fluorescence  by 
a  two-level  atom  in  a  monochromatic  laser  field.  There,  the  spectral 
distribution  of  che  fluorescence  consists  of  a  5-function  at  the  laser 
frequency,  superposed  on  a  smooth  background,  and  in  the  limit  where  equals 
che  laser  frequency  the  emission  rate  is  infinite.  In  any  case,  the  example 
of  this  section  shows  that  the  situation  of  perfect  resonance  has  to  be 
considered  with  caution. 


XV.  SCHRODINGER  PICTURE 


The  time  dependence  of  the  emission  rate  dW/dt  in  Eq.  (13.4) 
by  the  time  evolution  of  the  Heisenberg  operator  p( t) .  A  more 
representation  of  dW/dt  can  be  obtained  by  a  transformation  to  the 
picture.  When  we  cake  the  Schrbdinger  and  Heisenberg  picture  to 
t  -  0 ,  then  we  have 


is  governed 
transparent 
Schrbdinger 
coincide  at 


M(0)  -  M  -  JJ{  +  '>  +  ^ 


(15.1) 


with 
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M(+)  -  P  UP  ,  (15.2) 

-  g-  e 

M(*}  -P,P  .  (15.3) 

-  e-  g 

in  terms  of  the  projectors  P  and  P  onto  the  (possibly  degenerate)  excited 

eg 

and  ground  levels,  respectively.  In  free  evolution  (no  coupling  to  the 
radiation  field) ,  the  Heisenberg  representation  of  the  lowering  operator  p'  ^ 
is  f<^  +  ^(t)  -  exp(  -  i(jQt)u^  +  ^  .  which  equals  the  positive -frequency  part  of  p(t) , 
e.g.,  jj(t)^  +  ^  -  ^'  +  ^(t).  Due  to  the  interaction  with  the  electromagnetic 
field,  the  spectral  distribution  of  ^  +  ^(t)  will  acquire  a  finite  width  around 
the  central  frequencv  u^,  but  to  an  excellent  approximation  this  operator  will 
still  contain  positive  frequencies  only.  Hence  we  can  sec 


u(t) 


<±) 


(15.4) 


For  any  two  operators  A(t)  and  B(t)  we  have  the  identity 

<A( t)B( t)>  -  Trp( t)A(0)B(0)  ,  (15.5) 


in  terms  of  the  density  operator  p(c)  of  the  entire  atom  plus  field  system. 
In  Eq.  (13.4),  A(t)  and  B(t)  are  Cartesian  components  of  the  dipole  operator, 
and  in  the  Schrddinger  picture  they  act  on  wave  functions  in  the  atomic 
Hilbert  space  only.  Then  we  can  take  the  trace  over  the  field  states  in  Eq. 
(15.5)  according  to 


pA  t) 
a 


Tr 


fieldp(t) 


(15.6) 
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where  p  (t)  is  the  reduced  density  operator  for  the  atom,  irrespective  of  the 

cl 

state  of  the  field.  In  this  way,  we  find  for  the  emission  rate 


(-)  ..(+) 


(-)  .(+) 


Tr  p  (t)  lb|i£[j  •£;[  + 


+  c!|ifli+^  * ff n  ^  +  ^ 


(15.7) 


where  the  trace  runs  over  the  atomic  states  only.  Equation  (15.7)  expresses 
that  we  can  find  the  emission  rate  at  time  t,  once  the  atomic  density  operator 
_'_(t)  is  known,  since  the  operator  in  braces  is  completely  determined  by  Eqs . 
(15.2)  and  (15.3). 


XVI .  TWO- STATE  ATOM  AND  IDEAL  PC 

In  order  to  exhibit  the  principle  features  of  the  emission  of  phase- 
conjugated  fluorescence,  we  work  out  the  example  of  a  two-state  atom,  in 
combination  with  the  model  (ideal)  PC  from  Sec.  XIV. C.  A  model  two-state  atom 
has  an  excited  state  |e>  and  a  ground  state  |g>,  and  the  projectors  on  these 
states  are  Pg  -  |  eXe  j  and  -  |gxg|,  respectively.  Then  it  is  easy  to  work 
out  the  emission  rate  dW/dt  from  Eq.  (15.7),  and  we  obtain 


dW  1  0 

dc  “  2  A*“o(ne(C)  +  lpl  Vc>> 


(16.1) 


where 


a  -  t  l<e la  I &> I 

3re  p(cJ 


(16.2) 
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equals  the  Einstein  coefficient  for  spontaneous  emission  by  an  atom  in  empty 
space.  The  time  dependence  of  dy/dt  enters  through  the  populations  of  the 
atomic  levels,  defined  by 

n.(t)  *  Tr  p  (t)P.  ,  i  -  e,g  (16.3) 

i  a  a  i 

These  quantities  cannot  be  determined  with  the  present  theory,  since  the  time 
evolution  of  p  (t)  is  governed  by  the  coupling  of  the  atomic  dipole  to  the 
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vacuum  field  (and  possible  external  fields).  Also,  p  (t)  depends  on  the 
preparation  of  the  system  at  t  =>  0 .  Nevertheless,  when  for  a  given  time  t  the 

roouiations  are  n  t,t)  and  n  (t),  then  the  emission  rate  follows  from  Eq . 

e  g 

\16.1).  This  situation  is  analogous  to  the  fact  that  for  fluorescent  emission 
by  a  dipole  in  empty  space,  the  emission  rate  always  equals  dW/dt  -  Aj<wQne(t), 
no  matter  how  the  excited  state  became  populated. 

The  emission  rate  has  two  distinct  contributions.  The  first  term  is 
hAKw  n  (t),  which  equals  half  the  emission  rate  in  empty  space.  These  are  the 
photons  that  are  emitted  directly  by  the  atom  towards  tl  detector  (p-wave), 
and  the  factor  *s  accounts  for  the  fact  that  we  only  measure  the  emission  into 
the  half -space  z  >  0.  This  contribution  is  always  present  for  any  medium  in 
t  <  0,  simply  because  this  part  of  the  field  does  not  interact  with  the 
medium.  During  the  emission  the  atom  decays  from  its  excited  state  to  the 
ground  state  as  illustrated  in  Fig.  4. a.  The  other  half  of  the  number  of 
photons  which  are  emitted  in  the  |  e>  -  |  g>  transition  travel  in  the  negative 
z-direction,  and  they  serve  as  the  incident  field  on  the  medium.  For  linear 
reflection  this  would  give  rise  to  a  specular  wave,  with  intensity 
proportional  to  n  (t).  Interference  between  the  p-wave  and  the  r-wave  is 
incorporated  in  the  values  of  the  parameter  functions  bii  and  b^. 
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More  interesting  is  the  second  contribution  |  P  |  n^(t)  to  the 

2 

emission  rate.  This  term  is  proportional  to  j  P |  *”  and  represents  therefore  the 
phase-conjugated  fluorescence,  as  can  be  measured  by  a  detector.  The  question 
to  ask  is  what  physical  mechanism  is  responsible  for  the  emission  of  this 
radiation.  As  pointed  out  in  Sec.  XI,  an  interpretation  in  terms  of  an 
incident  wave  which  is  reflected  as  a  pc- wave  cannot  be  correct,  since  this  is 
in  conflict  with  the  different  retardation  times  of  the  various  photons. 
First,  we  notice  that  this  term  corresponds  to  the  actual  observation  of  a 
photon,  and  as  we  found  from  the  retardation  time,  this  photon  is  emitted  from 
the  site  of  the  atom.  This  emission  can  only  take  place  if  accompanied  by  an 
\e>  -  j g>  atomic  transition.  On  the  other  hand,  the  emission  rate  is 
proportional  to  n  (t)  so  that  the  atom  must  be  originally  in  its  ground  state. 
Consequently,  before  the  emission/decay  process  the  atom  must  be  excited,  and 
because  of  energy  conservation  this  can  only  happen  in  combination  with  photon 
absorption.  The  only  photon  source  in  the  problem,  assuming  no  external 
fields,  is  the  PC.  Therefore,  we  must  conclude  that  this  photon  is 
spontaneously  emitted  by  the  medium,  due  to  the  presence  of  the  dipole.  Such 
a  process  would  give  rise  to  the  emission  of  a  photon  with  frequency  u,  but  we 
recall  that  this  phase-conjugated  wave  must  have  frequency  2u>-u>o>  since  the 
incident  field  has  frequency  .  An  energy-conserving  process  then  requires 
that  a  second  photon  with  frequency  w  is  absorbed  by  the  atom,  which  gives 
rise  to  an  excitation  from  j g>  to  |e>.  Figure  4.b  illustrates  this  three- 
photon  process. 

XVII.  DEGENERATE  TWO -LEVEL  ATOM 

A  two- state  atom  description  is  an  idealization  which  is  commonly 
employed  to  study  the  principles  of  matter-radiation  interactions.  In  this 
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fashion  we  can  discover  the  fundamental  mechanisms  of  a  certain  process,  such 
as  illustrated  in  Fig.  4  for  the  present  problem,  although  the  model  is  not 
necessarily  verv  realistic.  In  this  section  we  consider  a  degenerate  two- 
level  atom,  where  the  excited  and  ground  levels  have  angular  momenta  jg  and 
j  ,  respectively.  These  levels  are  (2j  +1)-  and  ( 2 j  +l)-fold  degenerate,  and 
the  magnetic  states  are  indicated  by  |jeme>  anc*  I  j  gtDg>  °^v^ous  notation. 
Then  it  is  immediately  clear  that  a  two-state  description  cannot  account  for 
the  behavior  of  the  atom,  even  if  one  would  select  a  certain  transition 
between  two  states  bv  an  external  field.  For  instance,  in  a  j  ”  0>Je  ”  ^ 
system  we  can  select  the  m  =  «•  0  states  by  driving  the  atom  with  a 

1  i  near  lv  -  po  1  ar  i  z  e  d  laser  in  the  z-direction.  Stimulated  transitions  occur 
onlv  between  the  0  0>  and  | 1  0>  states,  and  in  addition  we  have  the  two 
spontaneous  processes  from  Fig.  4  between  these  states.  But  then  also  the  |l 
1>  and  |l  -1>  excited  states  have  a  nonzero  dipole  matrix  element  with  the 
ground  state  |0  0> ,  and  the  three-photon  loops  from  Fig.  4.b  will  generate 
fluorescence  whenever  the  atom  is  in  its  ground  state. 

In  this  section  we  evaluate  the  general  expression  for  the  emission  rate 
for  a  degenerate  two-level  atom.  in  terms  of  the  dipo  le  -  se  lec  t  ive  raising 
operator 


d  -  )  (1ralr|jm)|jra  XI  m  I 

r  l  Jg  g  IJe  e  J  e  e  Jg  g> 


r  -  -1,  0,  1  , 


(17.1) 


m  m 
e  g 


the  lowering  part  of  the  dipole  operator  attains  the  form 


(+)  ^  .t 


)  d;^  . 


(17.2) 
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in  terms  of  the  spherical  unit  vectors  e^  with  respect  to  the  z-axis.  The 

summation  over  r  —  ±  1  then  gives  fii  ^  ,  whereas  the  r  —  0  term  equals  the 

-i| 

perpendicular  part  of  +  ^  •  With  the  Einstein  coefficient  for  spontaneous 
emission  in  free  space  and  for  a  degenerate  system, 


A 


3  irt 


1  <  j  a  1  1  A*  I  |jg> 


and  with  the  change  in  notation 


b., 

I 


b 

T 


l  b 


X 


for  r  -  ±  1 


for  t  -  0 


(17.3) 


(17.4) 


C  i| 
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c  I 


for  t  -  ±  1 


for  r  -  0 


the  emission  rate  becomes 


dW 

dt 


"  )<w  ATr  p  (t) 

/  O  cl  cl 


(b  d  d1  +  c  d 1  d  ) 

T  T  T  TXT 


(17.5) 


(17.6) 


In  order  to  demonstrate  the  similarities  and  differences  with  a  two- 
state  atom,  we  consider  again  the  case  of  an  ideal  PC,  for  which  b  “  1  and  c^ 
-  |P|^.  With  the  identities 
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r  + 
)  d  d* 

L  T  T 


(17.7) 


T  *  2je  +  1 

L  drdr  ■  2j  +  1  Pg  ' 

r  s 


(17.8) 


in  terms  of  the  projectors  onto  the  two  levels,  we  find 


S  - 1  v 

g  6 


(17.9) 


'.-'ith  the  notation 


,  1  2je  *  1  |  ,2 

g  “  2  2j  +  1  A  |F| 
Jg 


(17.10) 


this  can  be  written  as 


^  An  ( t )  +  liu  A  n  ( t ) 
dt  2  o  e  o  g  g 


(17.11) 


where  the  populations  of  the  levels  are  again  given  by  Eq.  (16.3).  Comparison 
with  Eq.  (16.1)  shows  that  the  only  difference  is  that  the  ground-level 
emission  term  is  now  multiplied  by  the  geometrical  factor  (2j  +l)/(2j  +1). 
Notice  chat  only  the  total  populations  of  the  levels  determine  the  fluorescent 
emission  rate;  it  does  not  matter  how  the  populations  are  distributed  over  the 
various  states.  This  holds  only  for  an  ideal  PC,  and  not  in  the  more  general 
situation,  as  given  by  Eq,  (17.6). 


43 


XVIII.  CONCLUSIONS 

We  have  studied  the  emission  of  fluorescence  radiation  by  an  atom  in  the 

vicinity  of  a  phase - conj  ugating  surface.  With  a  decomposition  of  the 

particular  solution  (dipole  in  empty  space)  of  the  Maxwell -Heisenberg 

equations  into  an  angular  spectrum  of  plane  waves,  we  have  been  able  to 

construct  the  reflected  field  by  the  medium.  Under  the  assumption  that  the  PC 

operates  via  four-wave  mixing,  the  refit  ' 'd  field  could  be  expressed  in  terms 

of  the  classical  Fresnel  reflection  coefficients  R  and  P  .  These  parameters 

a  o 

incorporate  the  material  properties  of  the  medium  and  the  details  of  the 
geometry  for  the  four-wave  mixing  process  (like  the  polarization  of  the  pump 
beams,  and  the  layer  thickness).  This  procedure  yields  the  exact  solution  for 
the  electromagnetic  field  in  the  entire  region  z  >  0  in  terms  of  angular 
integrals,  and  with  an  asymptotic  expansion  we  have  derived  the  value  of  the 
field  in  the  far  zone,  in  terms  of  the  observation  angle  (&,<}>').  After 
integration  over  a  2 n  solid  angle  we  have  obtained  the  total  fluorescent  power 
in  the  half -space  z  >  0,  where  it  turns  out  that  this  power  can  be  expressed 
in  terms  of  four  geometrical  parameter  functions  b  j|  ,  bj_,  c  |  and  cj_,  and 
quantum  expectation  values  of  the  equal- time  autorcorrelation  functions  of  the 
Cartesian  components  of  the  lowering  part  of  the  dipole  operator.  By 
considering  various  limiting  cases  in  Sec.  XIV,  we  have  shown  that  the  terms 
proportional  to  b||  and  bj_  account  for  the  fluorescent  emission  directly  from 
the  atom  towards  the  detector.  Also  included  in  these  terms  is  the  specular 
reflection  of  photons,  and  it  appeared  that  both  waves  interfere.  This 
interference  is  again  a  pure  geometrical  and  classical  phenomenon,  and  is 
accounted  for  by  the  functions  b y  and  b^.  The  presence  of  the  phase- 
conjugated  wave  appears  as  a  separate  term  in  the  expression  for  the  emitted 
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power,  which  indicates  that  this  wav*  does  not  interfere  with  the  p-waves  and 
the  r- waves. 

In  working  out  the  model  case  of  a  two- state  atom,  we  have  been  able  to 
show  that  the  directly-emitted  photons  and  the  specular  photons  are  both 
produced  in  a  spontaneous  decay  of  the  atom.  Since  both  emissions  are  a 
result  of  the  same  decay  process  |  e>  -»  |  g> ,  these  photons  must  necessarily 
interfere,  as  is  expressed  by  their  joint  appearance  in  the  single  parameter 
functions  bn  and  b^.  On  the  other  hand,  the  probability  for  the  emission  of  a 
phase -conjugated  phocon  is  proportional  to  the  population  of  the  ground  state, 
which  renders  this  mechanism  independent  of  the  previous  one.  3v  considering 
the  retardation  times  of  the  various  photons,  we  can  track  down  the  underlying 
physical  mechanisms,  and  the  two  responsible  processes  are  illustrated  by  the 
diagrams  of  Fig.  4.  The  most  remarkable  result  is  that  an  atom  in  its  ground 
state,  and  close  to  a  PC,  can  effectively  emit  a  fluorescent  photon,  in 
addition  to  which  the  atom  excites  spontaneously.  We  have  interpreted  this 
phenomenon  as  a  result  of  a  three-photon  process,  as  depited  in  Fig.  4.b. 
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FIGURE  CAPTIONS 

Fig.  1 .  The  z  -  0  plane  separates  the  vacuum  z  >  0  from  the  nonlinear  medium 
in  z  <  0,  and  the  radiating  atom  is  located  in  x  -  y  -  0  and  z  -  h.  In  an 
angular- spectrum  representation,  the  dipole  radiation  equals  an  integral  over 
the  parallel  components  Ki|  of  wave  vectors.  For  z  >  h  and  0  <  z  <  h  these 
wave  vectors  are  K  and  K'  respectively,  and  their  directions  are  indicated  on 
the  left-hand  side  of  the  figure  for  the  case  of  travelling  waves.  The  jC- 
waves,  which  are  referred  to  as  p-waves  in  the  text,  travel  directly  from  the 
site  of  the  atom  in  z  -  h  towards  the  detector  in  z  »  h,  whereas  the  waves 
with  wave  vector  K'  serve  as  the  incident  field  on  the  medium.  In  the  case  of 
evanescent  waves,  the  field  amplitudes  decay  in  the  z-direction  as  shown 
pictorially  on  the  right-hand  side,  and  these  waves  travel  in  the  K| 
direction.  Also  shown  is  the  phase  convention  for  the  unit  polarization 
vectors,  both  for  K  and  K’  waves. 

Fig .  2 .  The  atomic  dipole  is  located  at  r  -  h,  and  the  emitted  fluorescence 
is  detected  within  an  angle  6  with  the  z-axis.  Both  the  directly-emitted  (p- 
wave)  and  phase -conjugated  (pc-wave)  radiation  appears  to  emanate  from  the 
location  of  the  dipole,  as  follows  from  the  retardation  time  r  -  (h/c)cos0  of 
the  p  and  pc-waves,  and  the  geometry  shown  in  this  figure. 

Fig,  3 .  The  retardation  time  in  the  observation  of  a  specular  wave  equals 
(2h/c)cos0,  as  compared  to  the  observation  of  a  p-  or  pc-wave.  The  geometry 
in  this  figure  shows  that  this  retardation  can  be  interpreted  as  the  delay 
time  for  a  photon  which  travels  first  from  the  atom  to  point  A,  where  it  is 
subsequently  reflected  according  to  Snell’s  law  in  classical  optics.  Also, 
the  distance  2hcos$  equals  the  separation  between  point  B  and  a  mirror  dipole 
m'  at  a  distance  h  below  the  surface.  For  a  detector  within  an  angle  8,  the 
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r-photons  appear  to  come  from  this  mirror  dipole,  as  can  be  seen  from  the 
geometry  shown  in  the  figure. 

Fie.  4.  This  figure  illustrates  the  two  atomic  processes  which  contribute  to 
the  observable  fluorescence.  Diagram  (a)  represents  an  atom  which  decays  from 
its  excited  state  |e>  to  its  ground  state  | g>  under  emission  of  a  photon  into 
the  direction  of  the  detector.  This  is  the  p-wave,  and  the  mechanism  is  not 
affected  by  the  presence  of  the  medium.  The  mechanism  for  the  emission  of 
phase -conjugated  fluorescence  is  shown  in  diagram  (b) .  The  atom  is  originally 
in  its  ground  state.  Spontaneous  emission  of  a  photon  (pc -wave)  by  the  PC, 
and  subsequent  absorption  of  this  photon  by  the  atom,  brings  the  atom  to  its 
excited  state.  Spontaneous  decay  is  then  accompanied  by  the  emission  of  the 
observable  photon.  A  second  absorption  of  a  photon  with  frequency  w  then 
completes  the  process,  thereby  leaving  the  atom  in  its  excited  state. 
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